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The study of crowd dynamics is interesting because of the various self-organization phe- 
nomena resulting from the interactions of many pedestrians, which may improve or ob- 
struct their flow. Besides formation of lanes of uniform walking direction and oscillations 
at bottlenecks at moderate densities, it was recently discovered that stop-and-go waves 
[D. Helbing et ai, Phys. Rev. Lett. 97, 168001 (2006)] and a phenomenon called "crowd 
turbulence" can occur at high pedestrian densities [D. Helbing et al., Phys. Rev. E 75, 
046109 (2007)]. Although the behavior of pedestrian crowds under extreme conditions 
is decisive for the safety of crowds during the access to or egress from mass events as 
well as for situations of emergency evacuation, there is still a lack of empirical studies 
of extreme crowding. Therefore, this paper discusses how one may study high-density 
conditions based on suitable video data. This is illustrated at the example of pilgrim 
flows entering the previous Jamarat Bridge in Mina, 5 kilometers from the Holy Mosque 
in Makkah, Saudi-Arabia. Our results reveal previously unexpected pattern formation 
phenomena and show that the average individual speed does not go to zero even at 
local densities of 10 persons per square meter. Since the maximum density and flow 
are different from measurements in other countries, this has implications for the capac- 
ity assessment and dimensioning of facilities for mass events. When conditions become 
congested, the flow drops significantly, which can cause stop-and-go waves and a fur- 
ther increase of the density until critical crowd conditions are reached. Then, "crowd 
turbulence" sets in, which may trigger crowd disasters. For this reason, it is important 
to operate pedestrian facilities sufficiently below their maximum capacity and to take 
measures to improve crowd safety, some of which are discussed in the end. 

Keywords: Pedestrian dynamics; crowd turbulence; video analysis; fundamental diagram 

1. Introduction 

It is well-known that driven many-particle systems often constitute complex sys- 
tems, in which different kinds of pattern formation phenomena are observed [3]. 
Depending on certain system parameters (the "order parameters"), one may find 
transitions from one state of collective behavior to a qualitatively different behavior 
[5]. Such transitions typically occur when a certain "critical threshold" is crossed. 
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It is interesting to study, whether such transitions can also be found in systems in- 
volving humans |121l36j. In this case, the role of particles is replaced by individuals, 
which follow different interaction rules. Nevertheless, it has been shown that many 
stylized facts of pedestrian crowds can be well understood by so-called "social force 
models" [14] and other modeling approaches such as cellular automata [U [23] . As 
examples we mention the segregration of different walking directions into lanes with 
a uniform direction of motion, or oscillations of the passing direction at bottlenecks 
[19] . Recently, we have observed two unexpected transitions in extremely dense 
pedestrian crowds |17j . While we found laminar flows at small and moderate densi- 
ties, there was a sudden onset of stop-and-go waves, which was later replaced by a 
phenomenon of highly irregular motion called "crowd turbulence" . This dynamics 
is dangerous, as it may cause people to fall, and it seems to be related with coordi- 
nation problems between neighboring pedestrians competing for little space. Both, 
stop-and-go waves and crowd turbulence were previously not expected, because the 
acceleration time of pedestrians is only about 0.5 seconds i28j, which suggests a 
quasi-adiabatic relaxation to a stationary state characterized by the "fundamental 
diagram" , i.e. a flow-density relationship. In contrast, stop-and-go waves in freeway 
traffic are caused by significant delays in the speed adjustments of vehicles [13) . 

The study of dense pedestrian crowds is particularly interesting, since it is one of 
few systems involving a large numbers of humans (up to millions) , where the detailed 
dynamics can still be revealed by analyzing empirical data [17j or performing agent- 
based simulations |i43j . Particularly intensive research activities have been triggered 
by the study of stampedes [16] . While "panic" has recently been studied in animal 
experiments with mice [34] and ants [3] there is still an evident lack of data on 
critical conditions in human crowds. Recent work f35', '4? goes into the direction 
of comparing empirical studies with each other, and bringing consensus to some of 
the fundamental issues in evacuation dynamics. 

However, for a long time, unidirectional pedestrian flows were predominantly as- 
sumed to move smoothly according to the " fluid-dynamic" flow-density relationship 



where Q represents the flow per meter width, p is the pedestrian density, and the 
average velocity V is believed to go to zero at some maximum density as in traffic 
jams [5j [33ll22ll32l im [37] . Formula ((U is often used as a basis for dimensioning and 
designing pedestrian facilities, for safety and evacuation studies. However, empirical 
measurements are often restricted to densities up to 4-6 persons per square meter 
only. In the study [41] . for example, the maximum density Pmax is 5.4 persons per 
square meter, and the corresponding flt curve of the speed-density relationship is 



where Vq — 1.34m/s is the free speed at low densities and a = 1.913 persons per 
square meter a fit parameter. Some other measurements of pedestrian densities. 



m 



Q{p) = pV{p) , 



(1) 
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however, reach upto 6 persons per square meter or more [33[ I27j . So, which mea- 
surement is correct? And what happens at even higher densities? Based on the 
projected area of human bodies, upto 11 persons may fit into one square meter [38] . 

In order to address these issues, we wih evaluate video recordings of the annual 
Mushm pilgrimage, where conditions are known to be particularly crowded. We will 
also see that these data are not directly transferable to Western European condi- 
tions and vice versa. Video-based studies of pedestrian flows have recently been 
carried out by other authors as well [26l [20l [37] . Most of them focus on bottleneck 
rather than unconstrained flows. Helbing et al. [H] have recently proposed a macro- 
scopic model for such flows. The situation we analyze in the following, however, is 
characterized by particularly extreme densities, which have finally led to the crowd 
disaster on January 12, 2006. Therefore, we expect particular insights into critical 
crowd conditions. 

Our paper intends to identify some reasons for the different flow-density curves 
(fundamental diagrams) reported in the literature and to close the data gap in the 
safety-relevant range of extreme densities. Considering the inconsistent measure- 
ments in the literature, which one should a capacity analysis be based on? Using 
the wrong curves can easily imply a wrong dimensioning and, thereby, future disas- 
ters. Generally, the difference between local and global measurements has not been 
sufficiently paid attention to. The same applies to the body size distribution and 
the cultural backgrounds (regardings prefered spacing and speeds). Many papers 
also don't present standard deviations, which are important to judge the variability 
of the measurement points, to compare different data sets and to determine the 
required capacity reserves when planning pedestrian facilities. 

The technique of video analysis is an important, but not the main point of 
our paper. Nevertheless, we show in detail how certain problems related to video- 
based evaluations can be successfully overcome and how the results depend on the 
specification of parameters. Note that the video-based evaluation method described 
in this paper was also used for another study [17]. While that study focussed on 
the dynamics of crowd disasters, this paper addresses the measurement process and 
safety-relevant features of the speed-density and fiow-density diagrams. These issues 
are important to draw correct conclusions from the video data. Another focus are 
the determination of critical crowd conditions. It will turn out that neither the 
density nor the speed or fiow field are good measures of the criticality in the crowd, 
as the latter actually depends on the crowd dynamics. 

The remainder of this paper descibes the large efforts and various difficulties 
that must be overcome in order to extract the characteristics and dynamics of 
crowd behavior at large densities from video recordings. As commercial tracking 
software did not do the job, we had to develop new algorithms which were capable 
of dealing with hundreds rather than dozens of people in a fully automated way. 
This involved the evaluation of terabytes of data, recorded at many measurement 
sites over several days in January 2006 (during the Hajj of the year 1426H). The 
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whole project, including the comparison with manual counts for validation and 
calibration of the method required several man years. 

Our paper is structured as follows: Section [2] will describe our technique of video 
evaluation, which is the basis of our data analysis. In Sec. [3] we will then define 
a new measurement method for local densities and flows. Problems like how to 
correct for the hiding of people by umbrellas will be addressed as well. Afterwards, 
Sec. E] will present empirical measurement results for the flow-density and speed 
density diagrams at the entrance of the previous Jamarat Bridge. Here, we will 
also study the considerable variability of density, velocity and flow data. In Sec. [5] 
we will describe the stop-and-go waves and the phenomenon of "crowd turbulence" 
discovered at high densities. Moreover, we will determine warning signs of critical 
crowd conditions, in particular the "crowd pressure" . Section [6] will summarize 
our results and discuss implications for the safety analysis and dimensioning of 
pedestrian facilities. 

2. Measurement Site and Video Tracking 

2.1. Description of Measurement Site and Conditions 

In order to get a better understanding of extremely dense conditions, we have 
scientifically evaluated the situation during the stoning ritual in Mina through 12 
fixed cameras mounted on high poles. The overall video material of the 10th to 12th 
day of Dhu al-Hijjah, 1426II amounts to more than 2 Terabyte of data. A special 
focus of the study presented here will be on the entrance area of the Jamarat Bridge 
(see Fig. [T]), where the situation became particularly crowded and the sad crowd 
disaster happened on January 12, 2006 (the 12th day of the pilgrimage). 



Figure 1. Illustration of the old Jamarat Bridge and the video-recorded area we are concentrating 
on in this study. There were many more cameras installed, but the one in the entrance area showed 
the highest densities and the most interesting crowd dynamics. 
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2.2. Video Tracking Method 

Video tracking has become a common and comfortable tool for empirical pedes- 
trian research recently [39l EH [25] . Its big advantage is the automatic evaluation of 
pedestrian trajectories, but video tracking has also its limitations: 

• Some tracking softwares require to specify the starting points of pedestrians 
manually by clicking on their heads. 

• Possible camera positions often imply small areas of recording or oblique 
camera positions. The latter implies that pedestrians may be hidden behind 
each other at high densities. In both cases, one must correct for effects of 
perspective. 

• In dense pedestrian crowds, pedestrians are often "lost" or interchanged by 
the tracking algorithm. 

• Most tracking algorithms are restricted to tracking several dozens of pedes- 
trians for the above reasons and reasons of numerical performance. 

Thanks to the 35 meter high pole at which our video camera was mounted, the 
recorded area was quite large (28mx23m) and the recordings could be made almost 
perpendicular to the ground. Therefore, pedestrians of different height were usually 
not hidden behind each other. However, the radius of a pedestrian head extended 
over 2-3 pixels only (corresponding to about 15-25 pixels covered by each head). 

The software developed by ourselves automatically determines heads by seaching 
for round structures. For this purpose, several digital filters (transformations) are 
successively applied to the video frames to enhance their contrast and identify the 
relevant structures (see Fig. [2]) . To gain a higher accuracy, the following additional 
steps are performed: 

• Double-checking of the identified circles by applying an Artificial Neural 
Network [21] trained to recognize heads. 

• Application of adaptive histogram equalization to compensate for variations 
in the light conditions. 

• Application of corrections derived from extensive manual count data (see 
Fig. Hand SeclLH]). 

Figure [3] illustrates how the video frames look like during several processing steps 
of the tracking algorithm. 

For each frame, the identification of heads yields the locations ri{t) of the pedes- 
trians i at time t. By comparing these with their local neighborhoods in the next 
video frame, one can determine their velocities Vi(t). The speed information is used 
to estimate the location of the pedestrians in the subsequent frame, which improves 
the accuracy of the tracking procedure. With a resolution of 25 pixels per meter 
and 8 frames per second, it is possible to determine even small average speeds by 
calculating the mean value over a large enough sample of individual speed measure- 
ments. 
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Figure 2. Illustration of how the video- analysis software operates. The raw video data are fed from 
the top, and each frame is passing a number of filters (actually more than shown over here), until 
the head locations and velocities are obtained with sufficient reliability in the end. 



2.2.1. Comparison of Automated and Manual Counts 

The tracking routine docs not only give good estimates of the local densities, speeds 
and flows. It is also suitable for counting people. Figure O compares automated and 
manual counts for one street leading to the Jamarat Plaza. Despite the difficulty to 
distinguish people walking in opposite directions, the reliability of the automated 
counts is reasonably high. The deviation is actually of the same order as the devi- 
ation between manual counts of two different persons. 



October 25, 2008 9:41 WSPC/INSTRUCTION FILE 

acs'fundanientarrevised2 



From Crowd Dynamics to Crowd Safety 7 




Figure 3. Illustration of video frames before and after certain transformations performed by our 
video analysis software. Upper left: Original frame. Upper right: Same frame, after the lense 
distortion was corrected for (see the straightened wall in the lower part.) Bottom left: The frame 
after edge detection and thresholding were applied. Bottom right: Frame superposed to the head- 
detection probability-density surface, with crosses on the most likely head locations as determined 
by an Artificial Neural Network (ANN) I24| . The probability density is determined by searching 
for head-like patterns in the current frame as well as extrapolating the locations of head detections 
made in previous frames, which results in a reinforcement when the same head is detected in many 
consequitive frames. A cross is displayed if the probability density is above a certain threshold and 
there is no more likely head location within a distance corresponding to one head diameter. The 
inaccuracies in the upper left corner are both due to the fact that the video is blurry in this region, 
but also because these pedestrians have just entered the video and therefore the reinforcement of 
the tracking algorithm has not yet taken effect. 
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Figure 4. Left: Part of a video frame used for manual counting and calibration. Right: Schematic 
illustration of the method used for the manual counting of pedestrians. The crowd video is played 
for 5 seconds in slow- motion at 1/lOth of the original speed, and two people are independently 
counting the number of pedestrians who cross a certain line within these 5 seconds. The counting 
is not done for the complete video recording, but for short video sequences separated by time 
intervals of 10 minutes each. 
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Figure 5. Flow and cumulative flow of a pedestrian road to the Jamarat plaza over 72 hours. Top: 
Comparison of automatically and manually determined values, using the methods illustrated in 
Figs. [2] and |4] Bottom: Same comparison, when more sophisticated classifiers based on Artificial 
Neural Networks |24l were used. 
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3. Measurement of Local Densities, Speeds, and Flows 
3.1. Data Evaluation Method 

From the locations r,;(f) of the pedestrians i at time t, we have determined the local 
density p{f,t) — pf{r) at a location r via the formula 

p{r,t) = pf (f) = -1_ ^exp[-||r,(t) - r\\^ / R^] , (3) 

where i? is a parameter. The local velocities have been defined via the weighted 
average 

Y.,v3^M-\\m-A\VR^ ] 

E,exp[-||rlW-f||Vi?2] 

and local speeds as v{f,t) = Vt{r) — \\vj^{r)\\, while local flows have been obtained 
according to the fluid-dynamic formula 

q{r,t)^q^{r)^ p{r,t)v{r,t). (5) 

The greater R, the greater the smoothing area around r. It can be calculated that 
the weight of neighboring pedestrians located within the area Aj^ = ttR^ of radius 
R is 63%. 

In principle, the average ("global") values g(t), V{t), and Q{t) of the density, 
speed, and flow, respectively, can be determined via the formulas for pf-{r), vf-{f) 
and q^{r), if the value of R is chosen sufficiently large. However, R must be signifi- 
cantly smaller than the radius of the video-recorded area in order to avoid boundary 
effects. Moreover, r must be chosen from the central part of the video so that the 
areas Aji around r are completely recorded. In practise, the global density is mea- 
sured by counting all pedestrians in the area of interest and then dividing by that 
area. The global velocity is determined by calculating the average velocity of all 
persons in that area, and the global flow is defined by the global density times the 
global velocity, or by the number of people passing a long cross section per unit 
time, divided by its length. 

In Fig. [6]we have studied by means of computer-generated distributions of pedes- 
trians, how well the average ("global") density values g are reproduced for different 
values of R. Our evaluation gives the following: The average {p^) of the local den- 
sity values pf according to Eq. ([3]) over different locations r agrees well with the 
average ("global") density g. Moreover, the variance of the local density measure- 
ments around the given, average density p, goes down with larger values of R. In 
fact, for R — > oo, all local density measurements result in the same value 

g lim pf[r,) . (6) 

R—^oo 

It corresponds exactly to the overall number of pedestrians, divided by the area 
Api = ttR^ they are distributed in, i.e. 
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Figure 6. Local density measurements according to formula (|3]l for 10 randomly picked points in a 
circle of 10 meter radius, compared to the average density in a circular area An of radius R = 100 
meters. Top: Uniformly distributed points with R = 1 (left) and R = 3 (right). Note that the local 
density of randomly distributed points varies strongly, as different points can be arbitrarily close 
to each other. Therefore, we also generated pedestrian distributions resulting from a pedestrian 
simulation with the social force model 1141 . which took into account the finite space requirements 
("diameters") of pedestrians. In this simulation, we assumed an average desired velocity of 1.34 
meters per second and a standard deviation of 0.26 meters per second as in Ref. 1411 . The desired 
direction of walking were assumed to be the same for all pedestrians. The initial distribution was 
randomly chosen, but the density measurement was made after a simulated time of 30 seconds. The 
resulting density distribution for R = 1 (bottom left) is significantly smaller, as the distribution 
of pedestrians is more regular than a random distribution due to their repulsive interactions. To 
a certain extent, it reminds of the distribution of points located on a hexagonal lattice (bottom 
right, for R = 1). 



Division of by Aji corresponds to the classical method of determining the aver- 
age density g, but it can also be obtained by averaging over local density measure- 
ments pf{r): 

i ' i 

(8) 
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Note, however, that we are not so much interested here in the average density 
but in the local density p{fi,t), as this is expected to determine the behaviour of 
pedestrian i. The variation of the values of Pi(r, t) in Fig. [6] is due to the statistical 
variation of pedestrians in space, which results in local density variations. If the 
value of R is fixed, the relative variation of measured density values is smaller for 
higher densities, due to the larger number of pedestrians in the area Nji = ttR^. 
This is a favorable property of our local density measurement method, as we are 
particularly interested in high densities. 

When measuring the local densities from video recordings of crowds, there may 
be additional reasons for density variations: 

(1) errors in the automated identification of pedestrians due to bad video quality 
(or, in cases of oblique cameras, due to hiding of other persons), 

(2) hiding of persons below umbrellas (to protect themselves against the sun), and 

(3) variations due to the spatio-temporal dynamics of the crowd. 

While the first problem can be partially healed by tracking over certain time periods, 
the second point will be addressed next. The third point was the subject of another 
publication [T7] . 

3.2. Dealing with Umbrellas 

Before we proceed with our evaluation, we have to discuss one particularity of our 
video recordings: A certain fraction p of pilgrims uses an umbrella for the protection 
from sun. These umbrellas mostly have a radius of 50 cm, so they can be used 
for calibration. However, they also cover an unknown number of pilgrims. In the 
determination of the average density g, one can correct for umbrellas as follows: 
Let A be the area covered by the video recording, Ai — A/N the area available for 
each of the N persons in this area, and A2 the typical area of an umbrella. Then, 
one umbrella corresponds to fc = A2/A1 covered persons. Moreover, let n be the 
number of visible people in the area not covered by umbrellas and m the number 
of umbrellas. The overall number of people is then N = n + km. If we assume that 
each p-th pedestrian carries an umbrella, we have m — pN, i.e. N — n + kpN or 
N — n/{l — hp) = n/(l — PA2/A1). The correct average density is g — N/A, while 
the density calculated without considering umbrellas is g' = n/A. Therefore, 



a.s g = 1/Ai = N/A. That is we would not have any corrections for p = 0, g = 0, 
or A2 = 0. However, in reality corrections are needed. These lead, on average, to 
higher density values. The corrected value for g can be calculated as solution of a 
quadratic equation. We get 



I-PA2/A1 l-pgA2 



(9) 



= 



1 - yi - ApA2g' 
2pA2 



(10) 
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While this solves the problem of how to correct the average density, if the fraction 
p and average area A2 of umbrellas is known, it does not help to correct the local 
density measurements, as the fraction of the measurement area covered by umbrel- 
las is varying strongly. This can lead to significant variations in the local density 
measurements, even if the local density is actually the same. 

For that reason, within a 20mxl5m area, we have randomly picked pedestrian 
locations and have determined the local densities around them. From this sample, 
we have removed the fraction 7 of lowest density values. Since the umbrellas are 
most likely to be found in the low-density regime, the probability of successfully 
filtering out umbrellas grows with 7. Note that we do not assume that persons with 
umbrellas tend to avoid large densites. Rather the idea is that the density calculation 
under- estimates the density in the vicinity of an umbrella because the surrounding 
persons are hidden under the umbrella. Fig. [7] shows for computer-generated data 
that filtering out 50% of the low-density data can still lead to significant under- 
estimation of the actual density by the measured one, while 7 = 0.95 leads to 
reasonably accurate density measurements. Therefore, in the following empirical 
evaluations of video-recorded pilgrim flows, we will restrict to sufficiently reliable 
measurements of local densities by using a cutoff value of 7 = 0.95. 




Figure 7. Simulated results to determine the influence of umbrellas on the local density measure- 
ments. The computer simulations have been carried out with the social force model, starting with 
random initial conditions. We have assumed that 2% of the persons carry umbrellas with a radius 
of 50 cm. While the actual local density has been determined from all generated points (repre- 
senting pedestrians), the measured local density has been determined by removing all points in 
a radius of 50 cm around 2% of the points. The figures show the measured over the actual local 
density for a cutoff of 7 = 0.5 (left) and 7 = 0.95 (right). A cutoff value of 7 = 0.95 reproduces 
local density measurements over a large range of densities well, while a cutoff of 7 = 0.5 tends to 
underestimate the actual local densities. 



Let us now check the plausibility of the above proposed methods to correct 
for umbrellas with real data. In order to investigate how the speed-density and 
flow-density diagrams depend on the cutoff value 7, we have evaluated them for 
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different values of 7 (see Fig.|8l). We find that, while the speed-density relationships 
vary relatively little with the variation of 7, there is a significant increase of the 
local flows for greater 7 values, particularly in the intermediate density range. This 
shows that fitting speed-density data may lead to unreliable conclusions regarding 
the flow-density relationship, while fitting flow-density data would lead to good 
velocity-density fits as well. The reason is that the density enters twice and in a 
multiplicative manner into the flow (as product of density and speed). 




02468 10 02468 10 



Local density (1/m^) Local density (1/m^) 

Figure 8. Empirical relationships between the local speed and the local density (left) and between 
the local flow and the local density (right) for the measurement site depicted in Fig.[l]and various 
values of 7 £ {0%, 25%, 50%, 75%, 95%}. The cutoff value 7 is supposed to correct for the influence 
of umbrellas. It is clear that the flow values are underestimated if pedestrians covered by umbrellas 
are not suitably accounted for, as for small values of 7. 

Moreover, we have determined the globally averaged speed and the global flow 
as a function of the global density for different assumed fractions of umbrellas p 
(see Fig.ll]). 

Let us now address the question, how the actual fraction of umbrellas can be 
determined from the video recordings. For this, we used the fact that the relationship 
between the local speed and the local density should not be changed significantly by 
the fraction p of umbrellas. We, therefore, defined a reference relationship pj-c{{v), 
measured the velocity v and density p, and calculated the fraction p of umbrellas 
that leads to consistent density values. In detail, the procedure was as follows: 

• As a reference relationship p^edv), we used a fit curve to data of the local 
density as a function of the local speed, determined for i? = 1 and 7 = 95%. 
According to our previous numerical studies, the large cutoff value 7 should 
eliminate the influence of umbrellas well. 

• At a given time point t, we picked 1000 random locations r^. 

• For these locations, we determined the local densities p'\^=o{fi) and the local 
speeds v\-y=a{fi) with R = Im and 7 = 0. The value 7 = ignored umbrellas 
and resulted in density values p' < p, but the speeds should be correct. 
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02468 10 02468 10 



Global density (1 /m^) Global density (1 /m^) 

Figure 9. Relationships between the globally averaged speed and the global density (left) and 
between the globally averaged flow and the global density (right) for various values of p G 
{0%, 1%, 2%, 6%}, using formula IIUI I. As a larger assumed fraction p of umbrellas implies 
a larger number of hidden pedestrians, it is clear that the flow must increase with the value of p. 



• We then estimated the corrected densities as p(ri) — Pi-cf{v\j=o{ri)) ■ 

• Next, we determined the global density Q{t) as average of the corrected local 
densities p{ri), i.e. g{t) = {p{fi)). Similarly, we obtained the global density g'{t) 
ignoring umbrellas as average of the densities p'|^=o('^i)i i-^. (}'{t) = (p'l^^ol'^i))- 

• Finally, we estimated the fraction of umbrellas via the formula 

which follows from Eq. ([9]). As the radius of most umbrellas was 0.5 meters, we 
used the value A2 = m? . 

The empirically determined fraction p{t) of umbrellas as function of time t is shown 
in Fig. [TOl It turns out that the fraction of umbrellas increases after 10:30am and 
reaches its maximum at noon time, when the sunshine is strongest. This shows the 
plausibility of our procedure, which has also been checked by manual counts. 



4. Empirical Findings 

After the previous description of our video evaluation procedure, let us now start our 
analysis of pedestrian recordings with a discussion of density measurements. Figure 
[TT]shows that there is always a large variation of local densities, and this variation is 
very important, as the safety in a crowd is not determined by the average density, 
but by the maximum occuring density. One can roughly say that the maximum 
densities are twice as high as the average density. Therefore, an average density of 
4 persons per square meter should not be exceeded [35]. 



October 25, 2008 9:41 WSPC/INSTRUCTION FILE 

acs'fundanientarrevised2 



Prom Crowd Dynamics to Crowd Safety 15 




% 10.5 11 11.5 12 12.5 
Time (h) 



Figure 10. Fraction p of pilgrims carrying umbrellas as a function of time. One can see that p is 
four times larger around noon time, as compared to an hour earlier. 




Figure 11. Distribution of local densities (with R = I) for a given average density (circles: 1.6 
persons/m"^, crosses: 3.0 persons/m^, dots: 5.0 persons/m^). Gamma distributions fit the his- 
tograms with 50 bins well (solid lines). After Ref. 1171 
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4.1. Relationships between Densities, Velocities, and Flows 

In this section we are mainly interested in local crowd conditions, as they are rel- 
evant for the criticality in the crowd. Let us first give empirical results for the 
dependence of the locally averaged speed v{r,t) = v^{r) and the locally averaged 
flow q{r,t) — q^{f) on the local density p{r,t) = p'^[r) (see Fig. [12]). We find 
that, with larger averaging radius i?, the maximum densities and flows are reduced. 
This is due to the considerable variation of the local density, speed and flow val- 
ues. In Fig. 1121 we have therefore also determined the density and flow values with 
a variable radius R = yJlQ/ q, in order to average over a comparable number of 
pedestrians in all density ranges. 




Figure 12. Average local speed (left) and flow (right) as a function of the local density p. The 
unexpectedly large range of local densities of up to 10 persons per square meter and higher are 
in agreement with manual counts obtained from photographs. Curves are plotted for constant 
R = Im, R = 2m, and for the density-dependent specification R = corresponding to a 

constant expected number of pedestrians in the area A^ = ttR^. One can see that smaller values 
of R tend to imply larger flow values, as the averaging is performed over smaller areas so that 
extreme values are not averaged out. In other words: The variation of local densities is smaller the 
larger the value of R. As cutoff value we have used 7 = 0.95. 

An important observation is the fact that the average local speed at the entrance 
to the Jamarat Bridge does not become zero even at local densities of 10 persons 
per square meter. This, of course, does not mean that pedestrians never stop, but 
that stops last for short time periods only (see Ref. [IT])- The observation of non- 
vanishing speeds at extreme local densities is in marked contrast to vehicle traffic, 
where drivers stop and keep enough distance to avoid collisions. The consequence of 
this pedestrian behavior is that their average flows remain finite, with short inter- 
ruptions only. However, there is a significant breakdown of the flow by a factor of 3, 
when the situation becomes congested. This breakdown implies a serious reduction 
of the effective capacity, which causes a further compression of the crowd until the 
situation becomes critical. In the worst case, this can lead to crowd accidents. 

In order to avoid a breakdown of pedestrian flows and over-crowding, the capac- 
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ity of a pedestrian facility should not be fully utilized. One rather needs sufficient 
capacity reserves to guarantee safety with respect to variations in the flows. From 
queuing theory it is known that huge queue lengths and enormous waiting times re- 
sult, if the inflow comes close to the flow capacity. However, when the waiting times 
become too long, people become impatient and pushy, which further deteriorates 
the situation. 

4.2. The Fundamental Diagram and its Comparison with Other 
Measurements 

Note that the actual capacity of a pedestrian facility can not be determined from 
Fig. 1121 as the maximum of local flows is much higher than the maximum of the 
average flows, i.e. the capacity is significantly below the maximum of this curve! 
The relationship between the average flow as a function of the average density is, 
therefore, displayed in Fig. 1131 This curve is called the "fundamental diagram". Its 
maximum flow is often used for capacity assessments. 

In order to compare our measurements with previous studies, we have digi- 
tized corresponding data published in the pedestrian-related literature. Figure [13] 
presents a comparison with various measurements. It turns out that the maximum 
average ("global") densities are higher and the average speeds are lower than the 
ones reported in many publications. Obviously, it is important to identify the rea- 
sons for this. A closer analysis shows that most data displayed in Fig. [13] are for 
European or North- American countries, while Mori and Tsukaguchi's measurement 
was carried out in an Asian country, where people are smaller. Therefore, the body 
size distribution has a dramatic influence on the velocity-density relationship and 
the fundamental diagram [31\ 1391 130j . It is consistent with observations that the av- 
erage "diameter" of pilgrims is smaller than that of European or American citizens, 
which explains the higher maximum densities observed in our study. 

It would be desireable, if measurements for different countries could be mathe- 
matically transformed to a universal curve, and if country-specific diagrams could 
be derived from it by determining a few parameters only. For example, one could 
think of using a density definition based on projected body areas, as proposed in 
Ref. 33J. This approach, however, is limited for the following reasons: 

(1) Densities still vary considerably, when the projected body areas have already 
reached 100% spatial coverage. 

(2) The coverage of two-dimensional space by projected body areas is hard to mea- 
sure, as it requires perpendicular recording, which is possible only in a small 
number of cases. 

(3) Favorable distances and speeds depend a lot on the cultural settings. For exam- 
ple, Ref. [12] shows that there are tremendous differences in walking speeds for 
different populations: The fastest average walking speed (Singapore) is three 
times higher than the slowest one (Malawi). During the Hajj, the low average 
speed at small densities is due to the fact that pilgrims tend to walk together 
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Figure 13. Global speed V (left) and global flow Q (right) as a function of the global density g in the 
video-recorded area. Our own data were determined using the time-dependent values p{t) depicted 
in Fig. llOl in order to correct for umbrellas. The average speeds and flows, obtained by averaging 
all data points for the same density, are represented by solid lines with error bars corresponding 
to one standard deviation (Note that most publications on the fundamental diagram of pedestrian 
flows do not provide error bars which must be criticized for neglecting the considerable variety of 
flow values that is also known for vehicular traffic). While variations of the average speeds appear 
to be relatively low, the flow values vary significantly, which must be taken into account by 30—40% 
safety margins in any capacity assessment. Symbols correspond to the empirical data of Mori and 
Tsukaguchi 1271 (circles), Polus et al. 1321 (squares), Fruin et al. [S] (triangles), and Seyfried et al. 
1371 (dots). The solid fit curve is from Weidmann 1411 . Note that the data by Mori and Tsukaguchi 
were not averaged over a large area and, therefore, rather represent local measurements. Note that 
the global flows determined from our measurement were bounded by the capacity of the stoning 
ritual at the pillars of the Jamarat Bridge, i.e. the maximum global flows at 3 to 6 persons per 
square meter could potentially be higher (as the local flow-density curves suggest). It should be 
stressed that the measurements here are global, i.e. averaged over the whole measurement area. 
Due to the large variability of local densities (see Fig. 1111 . this gives significantly lower densities 
and flows as compared to the similarly looking curves for local measurements presented in Fig. 1121 
Fig. 1 of Ref. [17], and Mori and Tsukaguchi ^27^. 



in groups, and a considerable number of group members is 50 years or older. 

Simple transformations can certainly not account for the varity of cultural factors, 
but a rough approximation may be made as follows: Determine the average speed 
at very low densities by measurement or extrapolation, and identify the value with 
Vb . Calculate the expected speed as a function of the density p = xpmax as 

V{xp^^^) = VoFix) (12) 

or the flow as a function of the density as 

Q(a;pmax) = Xpyna.yiVoF{x) , (13) 

where 



F{x) = 1 — cxp 



-0.35 ( - - 1 



(14) 



corresponds to Weidmann's curve As the maximum density Pmax is hardly mea- 
surable directly, it is treated as a fit parameter somewhere in the range between 4 
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and 12 persons per square meeter. and determined by minimizing the deviation of 
(|12p or (|13p from empirical measurements for the country or event of interest. Corre- 
spondingly scaled results are shown in Fig. 1141 The same procedure could be applied 
to improved standardized functions F{x), which take into accunt particularities at 
extreme densities. 




Scaled global density Scaled density 

Figure 14. Left: Comparison of scaled empirical measurements of the speed-density relationship 
with Weidmann's curve I I14I I (dashed line, see main text for details): The fit parameters for Mori's 
and Tsukaguchi's data 1271 (circles) are Vq = 1.40 m/s and Pmax = 9.00 persons/m-^, for Polus 
et al. 32 (squares) they are Vq = 1.25 m/s and pmax = 7.18 persons/m'^, for Fruin et al. [5] 
(triangles) they are Vb = 1-30 m/s and pmax = 6.60 persons/m^, for Seyfried et al. |37| (dots) 
they are Vq = 1-34 m/s and pmax = 5.55 persons/m-^, and for our data (thin solid line with 
error bars) they are Vq = 0.60 m/s and pmax = 10.79 persons/m^. The rescaling was done by 
first identifying Vb, either directly from the data or by extrapolating the speed data towards zero 
density. Then, a least-square fit was performed to find the value of pmax giving the best match of 
Eq. I|14| l. Finally, the densities were scaled by pmax, while the speeds were scaled by Vb- It can be 
seen that scaled speed-density data from different places in the world are reasonably compatible. 
Right: The compatibility of flow-density data is less obvious, particularly at extreme densities: 
The smooth curve was derived from Weidmann's speed-density relation II14I I. while the other solid 
curve corresponds to our scaled global flow as a function of the scaled global density. The dashed 
line (highlighted by circles) represents the scaled local flow (for 7 = 0.8) as a function of the 
scaled local density with Vb = 0.86 meters per second and pmax = 10.57 persons per square meter. 
Therefore, it must be stressed that flow-density relations should not be derived from speed-density 
fits, but separately fitted to empirical fiow data. 



5. Warning Signs of Critical Crowd Conditions 
5.1. Transition to stop-and-go waves 

Our study of pilgrim flows has also revealed dynamical phenomena. As one of the 
videos provided on the supplementary webpage j45) shows, at high crowd densities 
there was a sudden transition from laminar flows to stop-and-go waves upstream 
of the 44 meter wide entrance to the Jamarat Bridge [T7] . These propagated over 
distances of more than 30 meters (see Fig. [T5|) . The sudden transition was related 
to a significant drop of the flow, i.e. with the onset of congestion [T7]. Once the 
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stop-and-go waves set in, they persisted over more than 20 minutes. 

This phenomenon can be reproduced by a recent model based on two continuity 
equations, one for forward pedestrian motion and another one for backward gap 
propagation [18]. The model was derived from a "shell model" (see Fig. [15]) and 
describes very well the observed alternation between backward gap propagation 
and forward pedestrian motion. 

5.2. Transition to "crowd turbulence" 

After the occurence of the stop-and-go waves and the related breakdown of the 
pedestrian flow, the density reached even higher values and the video recordings 
showed a sudden transition from stop-and-go waves to irregular flows (see Fig. 1161) . 
These irregular flows were characterized by random, unintended displacements into 
all possible directions, which pushed people around (see Ref. [17j for sample trajec- 
tories). This "crowd turbulence" caused some individuals to stumble. As the people 
behind were moved by the crowd as well and could not stop, fallen individuals were 
trampled, if they did not get back on their feet quickly enough. Tragically, the area 
of trampled people grew more and more in the course of time, as the fallen pilgrims 
became obstacles for others [17]. The result was one of the biggest crowd disasters 
in the history of the pilgrimage. 

How can we understand this transition to irregular crowd motion? A closer look 
at video recordings of the crowd reveals that, at the time when the phenomenon of 
"crowd turbulence" occured, people were so densely packed that they were moved 
involuntarily by the crowd. Recently, a computer simulation of this situation was 
made [43] . It managed to reproduce the main observations with an extended social 
force model, assuming that people would try to gain space when the density becomes 
very high. 

5.3. Measuring Critical Crowd Conditions 

For a successful crowd management and control, it is important to know where and 
when critical situations are likely to occur, although it must be always kept in mind 
that not all hazards to the crowd can be reliably detected. Visual inspection of 
surveillance cameras is not very suited to identify critical crowd conditions: While 
the average density rarely exceeds values of 6 persons per square meter, the local 
densities can vary considerably (see Fig. [TT]) . Moreover, the contour plot of the 
density and its temporal evolution do not give a precise answer, where and when 
the situation becomes critical (see Figs. [17] a-|-b): At a certain time, the density is 
high almost everywhere in the recorded area. However, continuous alarms as given 
by warning systems in the past are not very useful, since security forces cannot 
be everywhere. Therefore, it is important to focus their attention and activity on 
specific areas. 

Also the analysis of the velocity field is not particularly suited to identify critical 
areas and times. Figures [T7]c-t-d show the negative divergence of the time-averaged 
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r(m) 

Figure 15. Top: Long-term photograph showing stop-and-go waves in a densely packed street. 
While stopped people depicted relatively sharp, people moving from right to left have a fuzzy 
appearance. Note that gaps propagate from left to right. Middle: Empirically observed stop-and- 
go waves in front of the entrance to the Jamarat Bridge on January 12, 2006 (after TT), where 
pilgrims moved from left to right. Light areas correspond to phases of motion, dark colors to 
stop phases. The "location" coordinate represents the distance to the beginning of the narrowing, 
i.e. to the cross section of reduced width. Bottom left: Illustration of the "shell model" (see 
Ref. 1181 ). in particular of situations where several pedestrians compete for the same gap, which 
causes coordination problems. Bottom right: Simulation results of the shell model. The observed 
stop-and-go waves result from the alternation of forward pedestrian motion and backward gap 
propagation. 
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Figure 16. Long-term photograph of the phenomenon of "crowd turbulence". In the fuzzy area 
of this picture, people are moved into all possible directions by physical forces building up in 
extremely dense crowds. 

velocity field times the density which, according to the continuity equation, describes 
the expected density increase along the trajectory of a pedestrian^ Moreover, Figs. 
[T7b+f show the magnitude of the curl of the time-averaged velocity field measuring 
the vorticity. Again, these figures give no indication of the time and location of the 
sad crowd accident. 

The decisive variable quantifying the hazard to the crowd is rather the variance 
of speeds, multiplied by the density. We call this quantity the "(crowd) pressure "El It 
allows one to identify critical locations (see Fig. fTTk ) and times (see Fig.lTTh). There 
are even advance warning signs of critical crowd conditions: The crowd accident on 
January 12, 2006 started about 10 minutes after turbulent crowd motion set in, 
which happened when the "pressure" exceeded the (i?-dependent) value of 0.02/s^ 
[T7] (see Fig. [T7h). Moreover, it occured more than 30 minutes after the average 
flow dropped below the (also i?-dependent) threshold of 0.8 pilgrims per meter and 
second, which can be identified as well by watching out for stop-and-go waves in 
accelerated surveillance videos (played fast-forward) [17j . Such advance warning 
signs of critical crowd conditions could be evaluated on-line by a video analysis 

''Prom the continuity equation dp/dt + ^(pv) = follows {dp/dt)/p = —\/v, where dp/dt = 
dp/dt + V ■ Vp represents the density change in time in the moving system, i.e. along an average 
trajectory. 

''This "gas-kinetic" definition of the pressure is to be distinguished from the mechanical pressure 
experienced in a crowd. However, a monotonoxasly increasing functional relationship between both 
kinds of pressure is likely, at least when averaging over force fluctuations. 
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Figure 17. Illustrations of (a, b) the density, (c, d) the negative of the divergence of the velocity 
field, (e, f ) the curl of the velocity field, and (g, h) the "crowd pressure" . The figures on the left 
show the quantities as functions of space, while the figures on the right show the same quantities 
as functions of time. A color version of these figures can be found on the supplementary web page 
1451 . The spatial plots are time-averages over the period from 11:45 to 12:30 on January 12, 2006. 
The arrows show the velocity field, obtained by averaging over the time interval from 11:45 to 
12:00. Brighter values correspond to higher values. The dashed ellipse indicates the area where the 
crowd accident on January 12, 2006, started. 
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system. This would allow one to gain time for corrective measures like flow control, 
re-routing of people, pressure relief strategies, or the separation of crowds into 
blocks to stop the propagation of Shockwaves. Such anticipative crowd control could 
certainly increase the level of safety during future mass events. 

6. Summary and Discussion 

While most previous measurements of the fundamental diagram for pedestrian flows 
have been restricted to densities upto 4 persons per square meter, we have mea- 
sured local densities of 10 persons per square meter and more on the 12th day of the 
Muslim pilgrimage in 1426H close to the entrance ramp of the Jamarat Bridge. The 
occurence of such high densities relates to the body size distribution of pilgrims and 
implies the possibility of unexpectedly high flow values. Even at densities up to 10 
persons per square meter, the average motion of the crowd is not entirely stopped. 
This can lead to over-critical compressions. However, the density is not the opti- 
mum indicator of critical crowd conditions. The "pressure" , defined as the density 
times the variance of velocities, provides better and more specific information about 
critical areas and times. 

Our results have been obtained with a newly developed, powerful video analysis 
method for pedestrians, which is well applicable even to dense pedestrian crowds. 
The automated evaluation delivers pedestrian counts in different directions of mo- 
tion and measurements of local densities, speeds, flows, and pressures. This infor- 
mation can be used for an improved surveillance of the crowd and the identification 
of critical crowd conditions. One warning sign is the breakdown of pedestrian flows 
under congested conditions, which causes stop-and-go waves and a further compres- 
sion of the crowd [17 . Another, very serious indicator of criticality is the occurence 
of high values of the "pressure" , which relates to the occurence of turbulent crowd 
motion [17j . This dangerous dynamics of the crowd, which may cause people to fall, 
can be well seen in accelerated videos. 

6.1. Critical Discussion of the Obtained Results 

Some of our findings, including the extreme densities of up to 10 persons per square 
meter and more, the finite average speeds even at those densities, as well as the 
stop-and-go flows and turbulent crowd motion, may be hard to believe. Let us, 
therefore, in the following address some questions regarding the reliability of our 
study: 

(1) Reliability of video tracking: The evaluation of dense crowds is certainly a 
difficult challenge, and the tracking of individual people is often not possible 
over long distances. However, this is also not needed to determine the densities, 
speeds and flowsEI Our method is based on sophisticated image enhancement 

'^If a head is identified for tfie first time at a certain position, we say tiiat pedestrian i has a 
trackabiUty level of 1, and generally, if it holds for n successive frames, the pedestrian has a 
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and pattern recognition methods, and it was independently evaluated, i.e. we 
did not have any influence on the manual counting procedure. When varying 
the time interval T (in minutes) between manual counts, the variance of the 
relative error £ between automated and manual counts was given by 



with V = 0.0015±0.0003 and a correlation coeflicient 7^2 = 0.72. Therefore, the 
variance of the relative error grows with the sampling time interval T approxi- 
mately according to a diffusion law (which represents a "random walk" of flow 
measurements over time). For T = 10 minutes, we have a standard deviation of 
VO.OIS = 12 percent, which was considered to be acceptable by the evaluators. 
Note that this deviation does not only reflect errors in the automated counting 
procedure, but also variations in the measured quantities over the time interval 
between successive manual counts. Therefore, the error of the automated counts 
is much smaller, and the proposed safety margin of 30% should be large enough 
to take into account both, temporal variations in the flow and measurement 
errors. 

Furthermore, note that the "error bars" shown in Fig. [13] primarily reflect the 
spatial and temporal variability of the flow, and not the measurement error. For 
vehicle traffic, for example, it is a well-known fact that flow measurements at 
medium and high densities vary strongly in time [13j . It is, therefore, likely that 
all original measurements of pedestrian flows have a large degree of variability, if 
the flow measurements for a given density are not averaged over. This variability 
is a matter of the time gap distribution |131 114j , which in turn is expected to 
depend on factors like walking individually or in groups, or the heterogeneity 
of the crowd with respect to culture and age. 

The accuracy of the global density measurements can be estimated from the 
distribution of local densities (see Fig. [TT]) . When the density in an area 
of 20m X 15m = 300m^ is determined from local density measurements with 
R = Im and An ~ 3.14m^, we have about 100 non-overlapping local den- 
sity measurements. Averaging over them reduces the standard deviation by 



VTOO = 10 as compared to the one in Fig. [TTl which results in a statistical 

error of less than 5% for the global density. Compared to the variation of the 
flow, this is negligible, so that it is justified to drop density-related error bars 
in Figs. [13] and [m 

(2) Extreme densities: Local densities of up to 10 persons per square meter or 
slightly more were not only determined from our video analysis. By evaluation 
of photographs and independently from us, they were obtained before. More- 
over, considerations based on projected body areas suggest that densities upto 

trackability level of n. With this approach, an automatically identified head will be more certainly 
a real head, the higher the trackability level is. ft is possible to require a minimum trackability 
level rimin to match manual counts well. An average trackability level of at least 4 turns out to be 
a good value for generating stable measurements. 



£^DT 



(15) 
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Figure 18. Snapshots taken in the video-recorded area for different densities. It can be seen that 
the recordings were made from high altitude with little distortion to correct for. The pictures 
show that the distribution of people is not at all homogeneous. Instead, the local densities vary 
considerably. One can also see some umbrellas used for protection from the sun. 

11 persons per square meters would be possible in principle j38], which is con- 
sistent with our results. Furthermore, Mori's and Tsukaguchi's data [57] come 
close to our measurements, when the scaling analysis of Fig. [14] is applied. The 
importance of the body size distribution can also be seen in Refs. [5T1 [551 [5U] . 
(3) Flow-density relationship: More important than the maximum possible den- 
sities (which should be avoided anyway) are the maximum global flow (as it 
determines the capacity of a pedestrian facility) , and the density at which the 
flow starts to drop (which gives rise to congestion and a further increase of the 
density). It should be stressed that an occurence of the maximum flow should 
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be avoided for two reasons: first, in order to prevent a breakdown of tlie flow 
and second, because the density at wliicli tlie maximum flow is readied is not 
comfortable to pedestrians anymore [5]. 

Note that both, our focal and global flow-density relations look different from 
the fundamental diagrams reported in the literature. In order to check the plau- 
sibility of our flndings, we have scaled our data by the maximum density pmax 
(which depends on the projected body area) and by the average desired velocity 
Vq (which is smaller for pilgrims than usual as they often walk in large groups 
including people aged 50 years and older). Our scaled speed-density data shown 
in Fig. I14h agree well with other measurements, which indicates their compati- 
bility, while we find deviations for the flow-density data at spatial occupancies 
of 0.6 and higher (see Fig. [Tib). These deviations can have different origins: 
In some cases, fundamental diagrams have been determined by fitting speed- 
density data and multiplying them with the density. The resulting fit curve will 
not necessarily represent flow-density data well, as the flow-density relationship 
is sensitive to deviations in the speed-density relationship. Furthermore, the 
data for the Hajj are expected to differ from measurements for other situations 
in which dense crowds do not move (e.g. on the platform of a subway station or 
in front of a stage). Therefore, the question is whether the main distinguishing 
feature, the second increase in the flow at extreme densities, is an artifact of 
our evaluation method or represents a real fact. 

(4) Second increase in the flow at extreme densities: This increase results 
from the measurement of finite average speeds even at extreme densities. In 
fact, it can be observed in the evaluated video recordings that (apart from 
intermediate stops reflecting stop-and-go motion) people kept moving at all 
observed densities, even after the accident occured and some areas were difficult 
to pass. People try to get out of areas of extreme density — this is why the 
average speed stays finite. 

It should be noted that the second increase of the fiow has important implica- 
tions for the dynamics [2]. In particular, it implies the coexistence of forward 
and backward moving shock waves. The minimum of the fiow-density curve at 
finite densities and crowd turbulence arc two closely related properties. And as 
crowd turbulence is well visible in the video recordings, when played in fast- 
forward mode (see the videos at http: //www. trafficforum.org/crowdturbulcnce), 
a second increase in the density is actually expected to exist. 

(5) Stop-and-go motion and crowd turbulence: Both phenomena have not 
only been observed once and in a single location. Stop-and-go waves can occur 
in large spatial areas and persist over many minutes, and we have seen the 
phenomena in recordings of different locations. The same applies to crowd tur- 
bulence. It has not only occured in January 2006, but also in previous years, 
and the observation is consistent with what Fruin has reported [6|: "At occu- 
pancies of about 7 persons per square meter the crowd becomes almost a fluid 
mass. Shock waves can be propagated through the mass, sufficient to ... propel 
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them distances of 3 meters or more. ... Access to those who fah is impossible." 
This visuahzes the conditions in extremely dense crowds quite well, and it is 
compatible with the location of the minimum in our flow-density data. 

(6) Local vs. global dependencies: It is clear that locally measured quantities 
vary more than global ones. If one averages over many local measurements, 
the maximum density and speed or flow are, of course, reduced. Nevertheless, 
according to Fig. 114b . the local and global curves agree surprisingly well, when 
the densities and speeds are scaled by the respective maximum values Vq and 
Pmax- As the local and global quantities have been determined in different ways, 
this is a good reason to trust the corresponding data. 

(7) Alarms of critical crowd conditions: The Hajj is in many ways an ex- 
treme event, with potentially critical situations occuring at different locations. 
If alarms are given whenever the global density exceeds a certain density defined 
to be critical, it is difficult to take specific actions. Such alarms may continue 
over many hours. The pressure-based method, in contrast, would give a warn- 
ing of potentially turbulent crowd motion over a short time period, and would 
indicate critical locations. It is no problem if the critical pressure (here: 0.02/s^) 
is shortly crossed before crowd turbulence takes over in a large area. This gives 
security forces a bit more time to prepare counter actions. 

Note, however, that counter actions should be taken much earlier, at latest when 
stop-and-go waves set in, as this indicates a breakdown of the flow. It should also 
be underlined that the critical flow of 0.8m/s [17] and the critical pressure of 
0.02/s^ depend on the parameters of the measurement method, in particular the 
choice of R, and it may depend on other factors like the kind of crowd as well. As 
a consequence, these critical thresholds need to be properly calibrated. In fact, 
rather than an automated alarm system we propose to view accelerated video 
sequences, which allows one to easily discover areas where a smooth flow is dis- 
turbed, where stop-and-go waves appear or where even crowd turbulence occurs. 
Furthermore, it would be helpful to enhance surveillance videos by overlaying 
additional information like the intensity of the crowd pressure, as illustrated by 
the video at http://www.trafficforum.org/crowdturbulence/pressure_video.mpg 

6.2. Some Measures to Improve Crowd Safety 

In the interest of crowd safety, densities higher than 3-4 persons per square meter in 
large crowds and particularly the onset of stop-and-go waves or crowd turbulence 
must be avoided. Therefore, a combination of the following measures is recom- 
mended: 

(1) Design: The infrastructure should be designed in a way that no bottlenecks 
or objects (e.g. luggage) will obstruct the flow. This can be quite challenging, 
if the usage patterns are changing. In particular, it must be avoided that the 
outflow capacity is smaller than the inflow capacity. Accumulations of large 
crowds should be avoided. If this is impossible, the outflow capacities must be 
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dimensioned such that the time to evacuate the area is much shorter than the 
time to fiU it. (It should not exceed a few minutes.) Furthermore, a system of 
emergency routes or, more generally, a "valve system" is required to be able to 
reduce pressures in certain areas of the system, where needed. 

(2) Operation: The infrastructure should be operated in a way that avoids coun- 
terflows or intersecting flows. Even merging flows may cause serious problems. 
Flows of people should be re-balanced if there is a large utilization in certain 
parts of the system while there is still available capacity in others. Obstacles 
(including people blocking the ways) should be removed from areas intended for 
moving. Some ways should be reserved for emergency operation and protected 
from public access. 

(3) Monitoring: Areas of accumulation or any possible conflict points (including 
crossing flows or bottlenecks, if these are really unavoidable in the system) 
must be monitored during highly frequented time periods. To support the job 
of the monitoring crew and maneuver security forces to the right places, it is 
helpful to display additional information in the surveillance videos, visualizing, 
for example, the density, the flow, and/or the crowd pressure. 

(4) Crowd Management: In certain events, the flow must be suitably limited to the 
safe capacity of the system (which should consider a safety margin of 30%). 
This may be done by applying a scheduling program, which is a plan regulating 
the timing and routing of groups of people. The compliance with the scheduling 
program must be carefully monitored (e.g. by control points and/or GPS track- 
ing), and deviations from it must be counter-acted (e.g. by fines). Moreover, 
an adaptive re-scheduling should be possible in order to respond properly to 
the actual conditions in the system. It is even more favorable to have a simu- 
lation tool for the prediction of the flows in the system. This would allow an 
anticipative crowd management. 

(5) Contingency Plans: For situations, where the system enters a critical state for 
whatever reason (e.g a fire, an accident, violent behavior, or bad weather con- 
ditions), one needs to have detailed contingency plans, which must be worked 
out and exercised in advance. 

The above is a list of only some of the measures that can be taken to improve crowd 
safety. For further measures see Refs. [38l[7l[29l[T0l[TTl|40l[9l[22]. 
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